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Lane-EmdenEquation: Picard vs Pade 
Zakir F. Seidov, Research Institute, College of Judea and Samaria, 

Ariel , 44837 Israel 

The approximate analytical solutions are found to Lane-Emden equation using 

Picard type iteration and rational Pade approximants . The presented values for 

the case of n=2 are more accurate than any one known in literature. 

Introduction 

The investigation of Lane - Emden (LEE) equation has a long - standing history 
(see, e.g., [1 - 10], to mention only part of a vast literature) . Recently 
Schaudt [6] has shown that Picard type iteration scheme may be used 
to show the existence, uniqueness and regularity of global solutions of LEE , 
in particular, for n > 1 . 

In this paper it was shown that straightforward usage of Picard iteration at 
general case of arbitrary n (even for n > 1) is not possible analytically, 
and so we apply this procedure for particular case of n = 2. 
The iterations are converging and allow to obtain the solution 
up to first zero of solution with great accuracy . 

Then we consider the modification of iteration scheme allowing 
to get very accurate Pade approximations and the values of relevant 
parameter at first zero of Lane - Emden function. 

We conclude that integral form of LEE allows to get (analytical) approximations 
in the terms of economized series (rather polynomials) and of Pade approximants 
in the way which is not straightforward from differential form of LEE . 

Presented results , for the case of n = 2, are more accurate 
than ones known in literature. 

ODE and integral equation 

Lane - Emden equation (LEEn) is a second order ODE with parameter n [1, 2] : 
Y" (x) + \y' (x) = -y n (x), y (0) = 1, y' (0) = 0. (1) 

By usual way (see, eg., [3-6]) it may be rewritten 
as nonlinear integral equation of Volterra type : 



y {x) = i 



s (i- -J) y n (s) ds. (2) 



This integral form of LEEn is more prof itable comparing to (1) in several aspects, 
reason for this being the explicit taking into account the initial conditions 
while deriving (2) from (1) . 

Iterative procedure in general case 

We use here (2) for obtaining approximate analytical solutions of LEEn . 
To this end, we rewrite (2) as : 
r x 



yjc+i (x) 



y k n (s) ds. (3) 



In this form , one clearly sees the possibility to solve the LEEn by Picard type 
iteration procedure , which should be converging at least for n> 1 [6] . 

Note that there is no such iterative procedure known for LEE in form (1) 

Now we start with y (s) - 1 . We get (for any n<.) : 
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yi (x) = l- — 



i 2 \ n 

■2 + n 3(2- -g-) i 2 I 3 5 



yz ( x ) = —■> + —n + "T x 2 Fi \-^, -n; — , —— 

1 + n (1 + n) 3 [2 2 6 

Next step leads to uncalculable integrals thus precluding any further 

(analytical!) iteration for general n, 

starting with ideal zero - step y„ (x) = 2 ( in the spirit of [ 6 ] ) . 

Case of n = 2 

So we consider particular values of n, e.g. n - 2 : 

2 4 6 

Y2 (x) = 1 -^ + ^o-T5ij- {6) 

Note that , in (6), only first three terms coincide with power series 
expansion of y (x) around x - 0. 

Rather each Picard iteration gives in some sense economized power series 
of the order which rapidly increases with iteration number. 

This allows to calculate function y (x) , and even the first zero 
of y (x) , with great accuracy . 

Second iteration, y 2 (x) (6), has the first zero X 2 - 3.89295, 
which is a rather rough estimation. In the Table 1, we present the values 
offirstzero, X ki such thaty^ (X k ) = , (for n - 2 ! ) 
from several Picard iterations : 

Table 1 . Parameters at first zero of solution to LEE2 

/ . X -> X pt 2 k-2 
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-x 2 y' (x) / 
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.89295362 


3.81298113844 
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.0983859052 
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4. 


.48911893 


4.40920514416 
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.2244983076 
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4. 


.33134231 


4.35230464399 
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.4143159741 


5 


4. 


.35410941 


4.35287035336 
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.4110793941 


6 


4. 


.35280006 


4.35281444601 


2 


.4110474718 


7 


4. 


.35281815 


4.35287459842 


2 


.4110459822 


8 




- 


4.35287463348 


2. 


4110465354 
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4.35287459594 


2. 


4110460122 


10 




- 


4.35287459595 


2. 


4110460121 


11 




- 


4.35287459595 


2. 


4110460121 



Picard iteration, power series and Pade approximation 

We observe the convergence of Picard type iteration scheme of solution to LEEn 
starting with zero approximation y (x) = 1, for LEEn of parameter n - 2. 
The reason of this convergence, one may guess, is that in Picard iterations, 
even in this rather simple case of LEE2, we get not exact power series, 
but in some sense economized power series . 

It is known [7], that exact power series expansion of function y (x) 
for arbitrary n has radius of convergence in all cases less than 
first zero of y (x) . So method of power series can not allow to find 
first zero of function y (x) . However one can use power series 
as an intermediate result to get from them Pade approximations, 
which are (much more) accurate than 
corresponding parent power series (see, e.g., [8, 9]). 

In order to illustrate the advantage of Pade approximation in comparison 
with Picard iterative scheme 

(not speaking of pure power series with their small radius of convergence ! ) 
we present so called diagonal Pade approximation of eighth order , PA[8, 8] , 
to solution of LEE2 : 



LEEPP.nb 



PA[8, 8} = 
(838857127878000 - 36744139139400 x 2 - 296471 637000 x 4 - 4578235200 x 6 - 39340559 x 8 ) I 
(6 (139809521313000 + 17177563695600 x 2 + 

483356654550 x 4 - 3069772860 x 6 - 18829469 x 8 ) ) . ( 7) 
It' s size is incomparably smaller than that of Picard iterations 
of the similar accuracy of evaluated value 
of first zero X Pr 8 = 4.35287035336. 

Pade approximation 

So we suggest the more effective modification of Picard iteration scheme. 
First we find power series using expansion of y " (x) 

in the integrand in (3) and then calculate the integral. We note that the 
general recursive relation for coefficients of power series for function 

y (x) of arbitrary n is given in [7] . Here we wish to underline that even 
not knowing these recurrence relations one can calculate the power 
series coefficients with iterative scheme from (3) while obtaining 

the coefficients of power series from (1) is not so straightforward . 

x 2 
We start with y 1 (x) = 1 - —— . 

6 

Then using iterative procedure (3) several times 

we get the power series for general n : 

x 2 n x 4 n (- 5 + 8 n) x 6 
¥s (X) = 1 ~ ~6~ + ~Y20 15120 + 

n (-3150 + 10805 n - 12642 n 2 + 5032 n 3 ) x 10 



3265920 1796256000 

n (138600 - 574850 n+ 915935 n 2 - 663166 n 3 + 183616 n 4 ) x 12 n x 14 



840647808000 1235752277760000 

(-21021000 + 101038350n - 199037015 n 2 + 200573786 n 3 - 103178392 n 4 + 21625216 n 5 ) + 

: " A (1891890000- 102 67435500 n + 23780949500 n 2 - 



1008373858652160000 



30057075285 n 3 + 21827357 636 n 4 - 8618115372 n 5 + 1442431856 n 6 ) . (8) 
Using this power series we get several Pade [k, k] approximations : 

ro „-, 60+ (-10 + 3 n) x 2 
PA[2 ' 2] = 60,3nJ > (5) 

PA [4, 4} = 
(45360 (-50 + 17 n) + 

420 (1250 - 951 n+ 178 n 2 ) x 2 + (-24500 + 29100 n- 10849 n 2 + 1290 n 3 ) x 4 ) I 
(45360 (-50 + 17 n) + 420 (350 - 645 n + 17 8 n 2 ) x 2 + 15 (190 n - 321 n 2 + 86 n 3 ) x 4 ) ; (10) 

PA[6, 6] {ao +S2 x2+a4 x4+ae X6) 



7 (b + b 2 x 2 +b 4 x 4 + b 6 x 6 ) 

a = 389188800 (12372500 - 28353900 n + 21372605 n 2 - 6679392 n 3 + 749132 n 4 ) , 
a 2 = 2494800 

(-592900000 + 1624109100 n - 1644080680 n 2 + 792866757 n 3 - 184792168 n 4 + 1 67821 56 n 5 ) , 
a 4 = 22680 (5780040000 - 18348690500 n + 23045300300 n 2 - 

14873877965 n 3 + 5246107880 n 4 - 963525716 n 5 + 72229136 n 6 ) , 
a 6 = (-3053351700000 + 10729693415000 n - 15644831 683500 n 2 + 12377539036650 n 3 - 

5750486943975 n 4 + 1567735653534 n 5 - 231 903413748 n 6 + 14339949064 n 7 ) , 
b = 55598400 (12372500 - 28353900 n + 21372 605 n 2 - 6679392 n 3 + 749132 n 4 ) , 
b 2 = 

3207600 (-30135000+ 98545300 n - 120932550 n 2 + 68800285 n 3 - 18368304 n 4 + 1864684 n 5 ) , 
b 4 = 3240 (807765000 - 3857983500 n + 7146037250 n 2 - 
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6578113455 n 3 + 3170490776 n 4 - 762978732 n 5 + 72229136 n 6 ) , 
b 6 = n (27533555000 - 124961404500 n + 221524322250 n 2 - 196938366435 n 3 + 
92676531282 n 4 - 21931005924 n 5 + 2048564152 n 6 ) . (11) 

Note that first two Pade approximations, PA[2, 2] andPA[4, 4], 

see (9) and (10) , are given first in [8] , 

where they are denoted ase^ 1 '' 1 (f) and©^ 2 ' 2 ' (£,) respectively, 

seeEq. (5) in [8] . The third Pade approximation PA [6, 6] , formula ( 11) , 

(and forth approximation given in Appendix) for LEEn solution of any 

parameter n are presented here at first. Also it can be shown that all 

1 ? 
four Pade approximations reduce to 1 - -—- x at n -> 0. Also at cases 

6 

n - 1 and n = 5 the formulas (9 - 11) and in Appendix 

provide Pade approximations of known functions : 

c-in I „\ / x 2 \ ~ 1/2 

(n = 1) , and \l + — - (n = 5) . 



x x " i 3 
In the last column of Table 1 the values of another parameter useful in the theory 
of polytropic stars (see [1, 2]) are also given according to Pade approximations. 

We note that values in Table 1 are of the largest accuracy known 
in literature (in [ 8 ] the values of X and \i are given with 7 decimals ) . 
The one important aspect of these calculations is using an ability 
of Mathematica [11] of exact calculations with integers, 
and obtained results in the last lines of Table 1 are even more 
exact and we omitted several last decimals in presented results. 
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Appendix 



PA\ 



(Co + C2 X + C4 X + Cs X + Cs X 



7 (do + d2 x 2 + d4 x 4 + dg x 6 + dg x 8 ) 

Co = 2778808032000 (- 405394 99623000000 + 

219612227381050000 n - 511477654026990000 n 2 + 676587517878560000 n 3 - 
562530934795015800 n 4 + 305353039511462775 n 5 - 108282312996133392 n 6 + 
24204097219027016 n 7 - 3097 60251290937 6 n 8 + 173153049930352 n s ') ; 

c 2 = 2724321600 (15975797005665000000 - 93812750271259250000 n + 

241429498430852550000 n 2 - 360556559054312350000 n 3 + 347275338859032654000 n 4 
225793301183157240975 n 5 + 100428697440662407020 n 6 - 30180690880628440984 n 7 + 
5866332008193485952 n 8 - 666264385302502832 n 9 + 33596198313370944 n 10 ) ; 

c 4 = 162162000 (-34824713815584000000 + 219927173906802600000 n- 
618139426187842370000 n 2 + 10251737437794 65588000 n 3 - 
1117966339422572119000 n 4 + 84297 9927142782418180 n 5 - 

448722785885998170189 n 6 + 168636577073120251356 n 7 - 43848080862147504184 n 8 + 
7512720295827146176 n 9 - 7 63484427 67 5133520 n 10 + 3487 4154904325696 n 11 ) ; 

c 6 = 83160 (3335453741881800000000 - 22401401379882945000000 n + 
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67833175458094929750000 n 2 - 122889080323677572450000 n 3 + 

148688095352232932540000 n 4 - 126727796509387793947000 n 5 + 

78043362960923937418225 n 6 - 34984902378314001772500 n 7 + 

11325758084851776233616 n 8 - 258129547 4254677262096 n 9 + 

393027792921365904112 n 10 - 35888044595615811200 n 11 + 1486082910296883968 n 12 ) / 
c s = -3935358993175533000000000+ 27489086733745765800000000 n - 

87587066649022538342500000 n 2 + 1 69052069768274 963800250000 n 3 - 
220942942891951760335450000 n 4 + 206673732164689256567505000 n 5 - 
142377428239163824382367250 n 6 + 7311126540259539569417 5125 n 7 - 
27964513897868903037647750 n 8 + 7861 679390615108424720096 n 9 - 
1577321939775971114089424 n 10 + 213635703656397601 638096 n 11 - 
17487522793907552014304 n 12 + 652628789624858161536 n 13 ; 
d = 2778808032000 (-405394 99623000000+ 219612227381050000 n- 

511477654026990000 n 2 + 67 6587 517878560000 n 3 - 

562530934795015800 n 4 + 30535303951146277 5 n 5 - 108282312996133392 n 6 + 

24204097219027016 n 7 - 3097 60251290937 6 n 8 + 173153049930352 n 9 ) ; 
d 2 = 8172964800 (3028027356585000000 - 18826223872160250000 n + 

51492765748754750000 n 2 - 81845560338319050000 n 3 + 83881693314626656000 n 4 - 

57961094822069523075 n 5 + 273402347 43773243460 n 6 - 86886647 84464616088 n 7 + 

1779913193666297344 n 8 - 212276122271447664 n 9 + 11198732771123648 n 10 ) ; 
d 4 = 32432400 (- 46946420101350000000 + 337879669634104500000 n - 

1084804099841582050000 n 2 + 2053550229663198700000 n 3 - 

2549884065663832883000 n 4 + 2182180740630633402950 n 5 - 

1313346140402699046975 n 6 + 555572535897326622972 n 7 - 1617 65775591138321896 n 8 

30859692537952223456 n 9 - 3470706639638745712 n 10 + 174370774521628480 n 11 ) ; 
d 6 = 83160 (283936435294050000000 - 2471215787875117500000 n + 

9596132853796098750000 n 2 - 22046318286704593700000 n 3 + 

33457928005411825895000 n 4 - 35407121047 497 988629250 n 5 + 

26826302480245895755525 n 6 - 14669178502968121829460 n 7 + 

5745112561862418989088 n 8 - 1571007 418300434014512 n 9 + 

284 607051855194 912368 n 10 - 30664360468491604352 n 11 + 1486082910296883968 n 12 ) ; 
d 8 = 9 n (17434185082476750000000 - 146647625747614792500000 n + 

551814890871351833250000 n 2 - 123195267 4494568458250000 n 3 + 

1822280070039019266490000 n 4 - 

1885364478216564017185250 n 5 + 14007704 99278411 609550225 n 6 - 

753286798160886430110090 n 7 + 290902568097394058143104 n 8 - 



78620673835912275662576 n 9 + 14104705308862137378704 n 10 
1507180671042570654496 n 11 + 72514309958317 573504 n 12 ) . 



